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ABSTRACT 

We provide estimates for the flux and maximum frequency of radiation produced when the 
magnetic field in a relativistic, highly magnetized, jet is dissipated and particles are acceler¬ 
ated using general considerations. We also provide limits on the jet Lorentz factor and magne¬ 
tization parameter from the observed flux. Furthermore, using the Lorentz invariance of scalar 
quantities produced with electromagnetic tensor, we provide constraints on particle accelera¬ 
tion, and general features of the emergent radiation. We And that the spectrum below the peak 
softens with decreasing frequency. This spectral feature might be one way of identifying a 
magnetic jet. 

Key words: radiation mechanisms; non-thermal - methods: analytical - gamma-rays: bursts, 
theory 


1 INTRODUCTION 

Relativistic jets where the energy is transported outward by Poynt¬ 
ing flux (ExB) have been invoked for many energetic outflows 
in astrophysical systems such as pulsars, quasars, micro-quasars 
and gamma-ray bursts (GRB). There is a vast peer reviewed litera¬ 
ture on this topic e.g., Michel (1969), Blandford & Znajek (1977) 
Blandford & Payne (1982), Kennel and Coroniti (1884), Begelman, 
et al. (1984), Coroniti (1990), Meszaros & Rees (1997), Lyubarsky 
& Kirk (2001), de Gouveia dal Pino & Lazarian (2005), Drenkhahn 
& Spruit (2002), Lovelace et al. (2002), Kulsrud (2005), Giannios 
& Spruit (2006), Komissarov et al. (2007), Tchekhovskoy et al. 
(2008), Metzger et al. (2011), Cerutti et al. (2012). 

The radiation is produced in these systems as a result of mag¬ 
netic field dissipation (referred to as reconnection, a generic phrase, 
which we will be using throughout this article), where particles are 
accelerated either via parallel electric field or first order Fermi pro¬ 
cess in converging flows, and they then emit photons via the syn¬ 
chrotron process. Radiation could also be produced in shocks inter¬ 
nal to the jet or when the jet interacts with the surrounding medium 
via a shock and transfers a fraction of its energy to particles in the 
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external mediunj^We do not consider the latter process in this pa¬ 
per. 

Magnetic reconnection is a complex and poorly understood 
process despite the work of numerous people on this problem for 
the last more than 50 years, e.g. Dungey (1953), Sweet (1958), 
Parker (1957), Petschek (1964), Syrovatskii (1981), Biskamp 
(1986), Yamada et al. (1997), Kulsrud (1998), Uzdensky & Kulsrud 
(2000), Bim et al.(2001), Drake (2006), Samtaney et al. (2009), 
Zweibel & Yamada (2009). Does this mean that we are doomed in 
our effort to understand those astrophysical systems where Poynt¬ 
ing jets play an important role until a predictive theory for recon¬ 
nections is developed? The answer depends on what it is that we 
want to understand about these systems. If we are interested in the 
general, global, properties then the fine details of the reconnection 
process might not matter. A basic understanding can be obtained 
from certain Lorentz invariant functions of electromagnetic tensor 
and conservation laws. The goal of this paper is thus modest, and 
highly restricted in this sense, i.e. to try to provide some constraints 
on Poynting jet parameters (without having to rely on a particu¬ 
lar reconnection model) so that magnetic dissipation can explain 


^ We are considering relativistic jets in this work which are Poynting flux 
dominated such as those that one encounters in Gamma-ray bursts, disrup¬ 
tion of a star by the tidal gravity of a massive black hole, or AGNs. If the jet 
energy were to be transported outward by particles as kinetic energy, then 
in that case the kinetic energy could be converted to radiation via internal 
and external shocks as discussed for GRB jets in the works of eg. Meszaros 
& Rees (1993), Rees & Meszaros (1994), Dermer et al. (1999), Ghisellini 
& Celotti (1999), Stem & Poutanen (2004), Beloborodov (2010), Thomp¬ 
son & Gill (2014); however the efficiency of converting jet kinetic energy to 
radiation in internal shocks is of order only a few percent eg. Kumar (1999). 
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some broad aspects of the data such as the efficiency for converting 
magnetic energy to radiation and the general shape of the emer¬ 
gent spectrum. In a recent paper Beniamini & Piran (2014) have 
provided constraints on a Poynting jet model for GRBs. Their gen¬ 
eral approach and results are very different from the one we pursue 
here. 

In section 2 we provide a few general properties of Poynting 
jet. We estimate the maximum energy electrons could achieve in 
reconnection, and the shape of emergent spectrum also in § 2 . 


2 POYNTING JET: A FEW GENERAL 
CONSIDERATIONS 

Figure [T] provides a schematic sketch of the system we are con¬ 
sidering. The magnetic fields of a relativistic Poynting jet undergo 
dissipation at some radiu^i? and a number of current sheets form 
within the causally connected region of the jet of comoving size 
R/V {R/V^ in lab frame) and efficiently convert the magnetic en¬ 
ergy to particle energy and radiation. Within any current sheet there 
are likely to be a number of different regions where particles are 
accelerated, and within a space of size R/T there are obviously 
many more. These acceleration regions are usually associated with 
X-points — located in between plasmoids or magnetic islands that 
form due to tearing instability - where the magnetic field vanishes 
in absence of a guide field and where the electric field can acceler¬ 
ate particles, or with converging flows where particle acceleration 
takes place via first order Fermi process. Regions where particles 
are accelerated will be referred to as PASs (particle acceleration 
sites). Some general considerations regarding particle acceleration 
in an individual PAS is discussed in §2.1| The maximum Lorentz 
factor (LF) of particles in PASs is determined by a combination 
of electric field strengtlj^and radiative losses in addition to energy 
conservation (O' While outside PASs, particles lose energy to 
radiation and any acceleration they experience is negligible. The 
particle distribution function inside PASs is not determined in this 
paper and is taken to be a hard powerlaw function as per numeri¬ 
cal simulations (Zenitani & Hoshino 2001; Jaroschek et al. 2004; 
Sironi & Spitkovsky 2011, 2014; Guo et al. 2014). However, the 
distribution function outside PASs is determined by solving an ap¬ 
propriate set of equations ( §2.2| l. Synchrotron radiation emanating 
from PASs and outside PASs is considered in §2.3[ In §2.4| we pro¬ 
vide an estimate for distance from the center where Poynting jet is 
likely dissipated and the number of current sheets in the causally 
connected region for an efficient conversion of magnetic energy to 
radiation. 

Let us consider a Poynting jet with magnetization parameter 
cr, Lorentz factor (LF) F, and isotropic equivalent luminosity L. 
The dissipation of magnetic field takes place when the jet is at ra¬ 
dius R, and that is also roughly the radius where radiation is pro¬ 
duced. The plasma is sufficiently cold before magnetic reconnec¬ 
tions so that the thermal pressure of particles can be ignored. 


^ The dissipation of Poynting jet could be spread over a wide range in 
radius — Ri to R 2 with R 2 ^ Ri . We are considering the maximum size 
of the region that is in causal contact at R 2 , i.e. between R 2 /C and R 2 . If 
this does not capture a good fraction of the energy dissipation process then 
we can add up results from other radii in a trivial way as processes going on 
in one region have no effect on another region that is not in causal contact. 
® Particles are also accelerated in converging velocity flows (first order 
Fermi process) and stochastic velocity fields (second order Fermi process), 
but these are not considered in this work. 




Figure 1. A schematic sketch of a Poynting jet and multiple reconnection 
zones within the causally connected region of comoving size R/V. Each re¬ 
connection zone has a number of particle-acceleration-sites (PASs), which 
are either regions between bright spots or plasmoids where magnetic field 
is small and particles are accelerated by electric fields or regions of con¬ 
verging flows where particles are accelerated by the Fermi process. The 
reconnection zones shown in this figure are artistic renditions of numerical 
simulations of Hesse & Zenitani (2007). 


The magnetic field in the jet comoving frame is Bq — all 
physical quantities in the jet comoving frame are denoted by a 
prime and observer frame variables are un-primed — which is re¬ 
lated to jet luminosity as 


L = B'o^T^R^c =» B'o= ={58G)^^^ ( 1 ) 

provided that a > 1 ; we are using the convenient notation = 
X/10". 

We adopt the standard model that charged particles are ac¬ 
celerated in reconnection layers where magnetic field dissipation 
takes place. The accelerated electrons with “thermal” LF 7 ' emit 
synchrotron photons of frequency less than or equal to v (in the 
observer frame) which is given by 


27rmec(l -I- z) ’ 


( 2 ) 


where q and me are electron charge and mass, and 2 ; is the redshift 
of the object. An upper limit to 7 ' can be obtained from energy 
conservation, i.e. the energy in accelerated particles cannot exceed 
the energy in magnetic field. This condition gives 


ne"/r, 


'■ {mp/me)a, 


(3) 


where is electron number density in the jet comoving frame, and 
cr = B'q / {S'Krif.mpc'^) (4) 


is jet magnetization parameter. The reason that equation 0 gives 
the maximum electron LF and not the average LF is because elec¬ 
trons accelerated in current sheets have a power-law distribution 
function {drie/d'ye oc 7 ^ 7 ^) with p < 2 , and therefore most of 
the electron “thermal” energy is carried by the highest energy elec¬ 
trons; numerical simulations find p < 2 when the region where 
the reconnection takes place is strongly magnetized, (t > a few, 
and when the reconnection layer is sufficiently large in size (e.g. 
Romanova & Lovelace 1992; Zenitani & Hoshino 2001; Sironi 
& Spitkovsky 2014; Bessho & Bhattacharjee 2012, Werner et al. 
2014, Guo et al. 2014). Making use of equations [T] & |^ for mag¬ 
netic field strength and electron LF, we obtain an expression for the 
maximum synchrotron frequency 


syn 

^max 


qK^^a'^{mp/me)^ 
2'KmeC^/'^R{l + z) 


= (2.2 X lO^eV) 


^ ^48 

.Ri5(1 + z) 


(5) 
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A more accurate estimate for that takes into account radia¬ 
tive losses are presented in §2.1| The synchrotron photons will be 
inverse-Compton (IC) scattered to higher energies by electrons pro¬ 
ducing these photons, and the maximum IC photon energy in ob¬ 
server frame is the smaller of nieC^fet ^/(I + ^) 3itd ~ vfmax- 
The specific flux at frequency v, i.e. flux per unit frequency, 
in the observer frame is 


/. 


cf B'^VNe \ z 
rrisC^ 4ndj ’ 


( 6 ) 


where iVe is the total number of electrons (isotropic equivalent) in 
the causally connected part of the jet with thermal LF ^ 7 e, and 
is the luminosity distance to the source. We can calculate the 
number of electrons needed to produce a given observed flux by 
combining equations 

iVe « 1.2 X 10 '^^ 28 f + 2)“\ (V) 


The optical depth of these electrons to Thomson scattering is, 
TT ~ = 8 X 10-^ + zr\ (8) 

and their “thermal” LF and kinetic energy luminosity they carry are 


2 'KmsvRc?^^ 

qLi/ 2 (i + 2 )-i 
iVe»TleC® 7 er 


1/2 


= 4X 10 " + 

jl/i ' 


{R/Vf 

= (1.2 X lO^^ergs-"). (lO) 


where UkeV is photon frequency (in units of 1 keV) for which the 
observed specific flux is fv,mjy Considering that the energy car¬ 
ried by electrons cannot exceed the energy in magnetic fields for a 
Poynting jet, we find 


Le/L<l 


r 7/12 


r<90 


(1 + 2 ) 


1/6 


f 1/3 - ,, 

Jv,mJyZ^18 


j1/6j2/3 


1/6 • 
feet7 


( 11 ) 


The reason for the approximate inequality sign in the above equa¬ 
tion is because magnetic fields of a Poynting jet could be com¬ 
pressed by a factor a few and thus Le could in principle exceed L 
by order unity. 

If we consider that there are rjp protons for every electrorj^that 
radiates at frequency v, then the kinetic energy luminosity carried 
by cold protons is 


Ner]pmpC^r 

{R/Tf 


(5 X 10^ ergs"^) 


— In I"' Vpf^^rnjyd]^ 28 


r 1/2 


( 1 + 2 ) 


.( 12 ) 


Therefore, the magnetization parameter for the jet at location where 
jet magnetic energy is dissipated and radiation is produced is given 
by 


^ ^ ^ 11 ** r-1 r3/2j-2 

If L.m.JyW8 rfL,28(l + Z)- (13) 

If 10% of electrons in the jet are accelerated, i.e. rjp = 10, and 
r = 20, then cr{R) ~ 25. And that means that the magneti¬ 
zation parameter at the let launching site where F ~ 1 is cro ~ 
r{R)a{R) ~ 500. 


^ 7 p > 1 when only a fraction of electrons in the jet are accelerated. 
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2.1 Particle acceleration in cnrrent sheets 


Consider an electron undergoing acceleration in a reconnection re¬ 
gion where the electric field is E', and the magnetic field is B'. 
In the absence of a guide field, the magnetic field vanishes at the 
X-point, and far away from it its magnitude is Bf but otherwise at 
this stage we place no further constraint on the electric and mag¬ 
netic fields. The electron starting from some place in the vicinity of 
the X-point is accelerated, and as it moves away it finds the strength 
of the magnetic field increasing. At some point when the magnetic 
field becomes sufficiently strong, which will be quantified shortly, 
the acceleration ceases if E' ■ B' = 0. However, even well be¬ 
fore this happens, the electron could stop accelerating due to ra¬ 
diative losses which will determine its terminal Lorentz factor. We 
consider this interplay between acceleration and radiative losses to 
determine maximum electron LF. 

It is best to view the motion of a particle acted upon by E' 
and B' from a frame where the fields point in the same direction 
(which is always possible except when \E'\ = \B'\ and the two 
fields are exactly perpendicular to each other). This special frame 
where E” || B” will be referred to as the AF frame (Aligned Fields 
frame). There are two quadratic Lorentz invariant functions of E' 
and B'-. 

h = -eci 3 ^sF°‘^E '^^/8 = E' ■ B' and 

I 2 = -F°‘^E^p/2 = E'"^ - B'^, (14) 


where Fa 0 is the electromagnetic tensor (anti-symmetric 2-form). 
Since E' ■ B' is Lorentz invariant, if there is a non-zero compo¬ 
nent of magnetic field in the direction of the electric field in one 
inertial frame, there will be a non-zero component in all inertial 
frames. However, the component of magnetic field perpendicular 
to the electric field can be made to vanish by frame transforma¬ 
tion if E'^ — B'^ > 0, or the electric field perpendicular to mag¬ 
netic field can be transformed away in an appropriate frame when 
E'^ — B'^ < 0. The point is that there exists an inertial frame (AF) 
where the transformed electric and magnetic fields are parallel, and 
the motion of the electron is the AF frame is as simple as can be — 
the electron momentum parallel to the fields increases linearly with 
time (if the electric field is non-zero in this frame), and the perpen¬ 
dicular component of momentum has a constant magnitude (time 
independent) and it rotates about the magnetic field at a constant 
rate. 

The simplest way to get to the AF frame is by a Lorentz boost 
in the direction E' x B'-, if E' x B' = 0, then obviously no Lorentz 
transformation is needed as we are already in a frame where the two 
fields are either parallel or one of them is zero. A straightforward 
Lorentz transformation algebra shows that the speed of the Lorentz 
boost required so that the fields are parallel in the new frame is 


Plt 


(1 -F e") - [(1 - e")" -F 4e^ cos" 6'] 
2e sin 9' 


where 


(15) 



cos o' 


E' ■ B' 

WW\' 


(16) 


Figure|^shows the LF of needed boost as a function of e for a few 
different values of 0'; a simple analytical expression for e ^ 1 and 
e < 1 is 


/3_lt ~ min |e, e "}sin6'', (17) 

which turns out to be exact (as opposed to approximate) for all 
values of e for the special case of 6 ' — -k 12 . 
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Figure 2. The Lorentz factor of the inertial frame (wrt jet rest-frame) in 
which the electric and magnetic fields in some region of current sheet 
are parallel to each other is shown in the upper panel as a function of 
e = E'!B'\ Tlt = (1 “ '^here Plt is given by equation 

{B} . The three different curves correspond to three different angles {O') be¬ 
tween electric and magnetic fields; O' = 70° (dashed cuiwe), 80° (dotted 
line), and 90° (solid line). Note that Tlj' ~ 1 except when electric field is 
almost exactly perpendicular to the magnetic field and the strengths of these 
fields are about the same. This makes it rather easy to can^ out calculations 
in the AF frame where E" || B" and transform variables back to the jet 
comoving frame. The lower panel shows the angle (measured in the jet co¬ 
moving frame in degrees) by which the electric field direction is rotated in 
the AF frame for three different values of O' which are same as in the upper 
panel; when the electric field vanishes is the AF frame then the angle is the 
rotation for magnetic field plus 7r/2. 


The electric field in the new frame follows from the two 
Lorentz invariant quantities mentioned above and is given by 


E”^ 


h + x/472 + q 
2 


and the magnetic field is 


(18) 


B" 


h_ 

E"’ 


(19) 


where Ii and I 2 are defined in equation l |14| l. The angle between 
the aligned electro-magnetic field in the AF frame and the electric 
field in the jet frame can be easily calculated and is 


cosO'p = 


(e — /3lt sin 0') 


( 20 ) 


E — 2ePLT sin 9' 

The lower panel of Figure]^ shows 9'^ as a function of t for a few 
different values of 9'. 

With these results in hand, we are ready to describe particle 


acceleration in a current sheet. Consider a charged particle in the 
vicinity of the X-point where E' B'. We can transform away the 
perpendicular component of the magnetic field by going to the AF 
frame, and in this frame the particle Lorentz factor 7 ” (the double 
prime emphasizes that we are in a different inertial frame, and not 
the jet comoving frame) increases as qE''t"/{msC?), which can be 
rewritten from the point of view of the jet frame as 


/ 

7e 


qeoB'ol' 

meC^ ’ 


( 21 ) 


where is the distance the electron has traveled along the direction 
of the electric field from its starting position in the jet comoving 
frame, and 


eo = E' jB'^. 


( 22 ) 


As the electron travels further and further away from the X- 
point, it feels the strength of the magnetic field increase and at some 
point when B' becomes stronger than E' the electron is no longer 
accelerated (unless E' ■ B' ^ 0) and its momentum vector gyrates 
about the magnetic field and the LF oscillations and drifts slowly 
with time. This generic behavior can been in figure 0 where nu¬ 
merical result for particle motion in a current sheet is presented. 

If the length of the region where E' > B' is £'^, then the maxi¬ 
mum LF of electron ~ geoBo£s/(meC^); it should be noted 

that (!^ 25 eo for a magnetic field configuration where B' increases 
linearly with distance from the X-point, and thus 7maa; S5 eo (see 
e.g. Larrabee et al, 2003). However, two effects can substantially 
limit electron LF below tbis value. One of which is “global" en¬ 
ergy conservation, which provides a limit for ^max as described by 
equation 0- And the other is radiative losses — synchrotron and 
inverse-Compton (IC) for systems of interest to us — that could 
restrict particle LF further. This is discussed below. 

Viewed from the AF frame where E” || B", the electron 
suffers radiative losses due to acceleration along the electric field 
direction, gyration about the magnetic field, and inverse-Compton 
scatterings. We evaluate each of these to determine the dominant 
loss mechanism, and its effect on 'y'^ax- The energy loss rate is cal¬ 
culated by first assuming that the magnetic field lines are parallel, 
i.e. B' ■ VB' = 0, and the electric field is nearly uniform. This 
estimate is then improved by relaxing these assumptions and by 
considering the more realistic possibility that magnetic field lines 
have non-zero curvature, and that the electric field has spacial fluc¬ 
tuations in the acceleration region. 

The power radiated due to particle acceleration along the elec¬ 
tric field can be calculated using the Larmor’s formula. The mo¬ 
mentum vector of the electron in the AF frame is nearly parallel 
to E” since it is being accelerated along the electric field and the 
magnetic field is parallel to E" in this frame. Therefore, the elec¬ 
tric field in the instantaneous rest-frame of electron is also E", and 
the magnitude of its acceleration in this frame is qE"/rris. It then 
follows from Larmor’s formula that the power radiated (which is a 
Lorentz invariant quantity) is otE"^ c/A -k ~ ut^B'^ cjA-n', from 
equation 10 and Fig. 1^ we know that F^t ~ 1 for e > 2, and 
hence E" ~ E' — soEq. This rate of loss of energy is indepen¬ 
dent of electron LF, and so the maximum electron LF in this case 
is bounded only by the size of the acceleration region. 

The synchrotron loss rate due to electron gyration about the 
magnetic field is crTB”^7”^/3x^c/(67r); where meC7”/3” is the 
4-momentum of the electron perpendicular to the magnetic field. 
In the region where particles are undergoing acceleration, the mag¬ 
netic field in the AF frame (B") either vanishes (if E'-B' — 0) or is 
parallel to E', and in either case the value of 7 ”/?” does not change 
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Figure 3. Acceleration of an electron in a cuiTent sheet is shown as a func¬ 
tion oftime in unit of l/cjc, where lJc = qB'^/{raec) is Larmor frequency. 
The electric and magnetic field configurations in the current sheet are taken 
from Larrabee et al (2003), i.e. sheet lies in the x-y plane with the electric 
field pointing in the x-direction and has a constant magnitude E' = cq-Bq, 
and the vector potential is A = — z^)/2tg]z\ we took eo = 0.01 

and the length of the “sheet” = 5 x 10^ for this calculation. The 

particle started out in the 2 = 0 plane with initial velocity of zero. The LF 
of the electron ( 7 ^) increases linearly with time as long as it is in the region 
where E' > B' (which is for about 2xl0^a;^^ for the parameters we have 
chosen for this calculation), and afterward when B' > E' the acceleration 
ceases and the electron gyrates about the magnetic field. These results are 
entirely consistent with analytical calculations presented in this section. 


ton energy density in the comoving jet frame i^t/j^Bo^/STr, and 
therefore the IC loss rate is 


dnieC^'y'^ 

dt' 


= fjTt/jsBo 7^^c/67r. 


(24) 


Equating the rate of energy gain for an electron as it is acceler¬ 
ated along the electric field with the rate of radiative loses we arrive 
at the following equation for the maximum value for LF 


+ iPb) 


Gtt 


qeoB'oC 

Or 

fmax ~ 1.5 X 10'^£^^(1 + 1p_ 


-l/2pl/2,-l/4pl/2 

B ] i 2 .^48 ^15 5 


(25) 


(26) 


where we made use of equation |Tto substitute for Bq, and we are 
considering the case where E' ■ B' 0. The maximum electron 
LF is the smaller of values given in equations 0 and \26\ . 

The distance an electron travels to get accelerated to 'ymax is 


^n. 


^meC 


geoBo 

r3/2,-3/4 3/2 

(4.7 X 10®cm) 7 , 

_ (5.2 X lO'^cm) 

The synchrotron photon energy corresponding to fru 

(150 MeV)reo(l + 

(200 eV) 


(27) 


(28) 


The minimum electron LF can be obtained by taking the 
length of the acceleration region to be no less than proton Larmor 
radius. This gives frnin ~ eo{mpfme). 


with time even as the electron continues to accelerate. Thus, the 
synchrotron loss rate (like the loss rate due to acceleration along the 
electric field) is nearly independent of electron momentum, which 
continues to increase linearly with time along E' while the electron 
is in the acceleration region. 

For realistic astrophysical systems we don’t expect the mag¬ 
netic and electric field lines to be perfectly straight in the accelera¬ 
tion region. The curvature of field lines and the variation of E'/B' 
with distance from the X-point causes the direction of E" to change 
(see Fig. I^for the dependence of 6 '^ on E' jB'), and therefore par¬ 
ticle momentum vector is also rotated. Due to these effects fj is 
no longer independent of time, and in fact even a modest curvature 
in E' would lead to fj ~ In this case the synchrotron loss es¬ 
timated above increases by a factor 7 '^ (the loss due to acceleration 
along E' increases by a similar factor) and is given by 

~ rrT(eo + sin^6l3)Bo^7'^c/67r ~ arB'^ cl&'K{2'i) 

where Bq sin 6 'g is the strength of the guide field. From here on 
we assume that the guide field is not much smaller than Bq and 
therefore particle acceleration is dominated by electric field parallel 
to the magnetic field. 

The inverse-Compton loss rate is proportional to the energy 
density of photons, which is closely related to magnetic field dissi¬ 
pation. Photons are produced via the synchrotron process in accel¬ 
eration regions and also outside it. Assuming that a fraction /pB of 
magnetic field energy in a causally connected region of size R/T is 
dissipated in a dynamical time and converted to radiation, the pho- 


2.2 Electron distribution function 


Simulations of particle acceleration in a reconnection layer show 
that the energy distribution is a hard powerlaw function below 
frnax exponentially cut-off above it, i.e. dn'f,/df oc 7 e~^°. 
for 7 ' < withpo < 2 . 

When we add up particle distribution functions in all PASs 
within the causally connected region of the jet at the radius where 
a good fraction of magnetic energy in the jet is dissipated, the re¬ 
sulting distribution is 


dn'f, 

dYe 


oc 7e 


for Ymin^fe^fp- 


(29) 


The value of p depends on how many electrons pass through PASs 
which can accelerate them to LF 7 e; if the number of PASs in¬ 
creases rapidly with decreasing j^ax then p > po. The electron 
distribution below fmin is either cutoff or drops off such that the 
total number of electrons with 7 ^ < 'y'min is small and can be ig¬ 
nored. The distribution above 7 ), also falls off more rapidly than 

77^. 

Let us assume that electrons spend an average of tY time in an 
acceleration region which is larger than I'a/c, the average is taken 


® Photons produced by the dissipation of magnetic field in a region of size 
R/r in the comoving frame travel a distance in a dynamical time which 
is also R/r, and hence all the radiative energy is confined to a volume 

~ _R3/r3. 
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over all particle acceleration sites or PASs in causally connected 
part of the jet. Furthermore, the total number of electrons injected 
into these acceleration regions, in the causally connected part of 
the jet, per unit time is A^e.cs- The average rate at which particles 
exit PASs should also he A^e.cs- Particles outside PASs cool down 
radiatively and therefore the particle distrihution outside is much 
steeper. We calculate this distrihution, and estimate the synchrotron 
flux from electrons inside and outside PASs. 

If the average time spent hy electrons outside PASs is then 
in that time electrons cool down to LF 


hyrmeC 


GTrmeC^F^i?^ 

(T’J'L 




-1 
48 ? 


(30) 


where 




R 


cTt' 


(31) 


is the ratio of dynamical time in jet comoving frame and . 

The electron distrihution outside PASs is obtained by solving 


djdK^^Jd-yj) d^UdK,dJdi,) 

dt' dj'e 

where the source function is 




(32) 


s'iO 


{p-l)Nl^cs f 7 . 


for 7min ^ 7e ^ 7p (33) 


p > 1, is the rate at which electrons with Lorentz factors ^ 
'ylnin leave acceleration regions and enter the surrounding medium, 
and 


_ 0/2 /2 
<^tBo 7e 

GnirieC 


(34) 


A quasi-steady state solution is reasonable to consider when the 
time it takes for a typical PAS in the jet to form and disappears is 
much shorter than the dynamical time, and there are many PASs 
in the causally connected region of the jet that contribute to par¬ 
ticle acceleration and radiation; the average of all these PASs can 
be taken to be roughly constant for about a dynamical time. The 
solution of equation l |32[ ( for p 1, in quasi-steady state, is easy to 
obtain and for 7 ' < 7 ^*^ ** given by 


rf7e 



'ymin ^ 7e ^ 7p 

(35) 

7c ^ 7e ^ ^min 


The above derivation assumes 7 ^ ^ 7 ', which should be a good 
approximation considering that 7 ^ ~ ^max ~ (eq. [ 2 ^ and 
7 c ~ 1 (eq. 30 1 . The distribution is effectively cutoff above 7 ^ 


outside the PASs since electrons of this high energy cool efficiently 
and their LF drops below 7 ^ quickly. The distribution function for 
the case where 7 p > 7 c > 7 min is 


dN' 


rf7e 


zKcs 


7c ^ 7e ^ 7p 


'ymin ^ 7e ^ 7c 


(36) 


Forp < 1, the source function is 
(i-p)iv;cs Ac 


s'W.) 


7p 


7p 


for ^ 7A 7p (37) 


and therefore most of the electrons are at 7 ' ~ 7 p. The solution of 
equation l|32f using the above source function for 7 ' < 7min is 


dK,d. ^ AA.cc ' 

die ^ 7c V7c/ 


for ie<ie^ip- (38) 


And the distribution function when ip > i > imin is given by 


die 



7c < 7 A ip 

(39) 

^min ^ 7e 7c 


The apparent discontinuity of the distribution function in equation 
l |39| > at 7 e = 7 c is because the two branches of solutions are inac¬ 
curate as 7 c approaches i. However, a steep drop off of the dis¬ 
tribution function just below 7 ' is physical. This is due to the fact 
that electrons with LF ~ 7 ^ (which are a majority of the electrons 
entering the medium in between PASs when p < 1) radiatively 
cool down to LF ~ 7 ' in the available time t'^^, and hence there 
is an accumulation of electrons in the neighborhood of i and that 
is responsible for a drop in the distribution function just below this 
LF. 


2.3 Synchrotron and IC spectra 


Electrons inside PASs gain energy due to acceleration along electric 
fields or as a result of Fermi mechanism operating in a converging 
flow field. The balance between energy loss and gain determines 
the terminal Lorentz factor for particles. Moreover, the particle dis¬ 
tribution function, and index p, are also determined by the accelera¬ 
tion and radiative loss processes, and thus the synchrotron spectrum 
due to radiation from electrons inside PASs is fi, oc for 

v'Si.i’p’, where 


qB'oi^V ~ 

27rmec(l -I- z) 27rmeC^/^i?(l -I- z) 


(40) 


is synchrotron frequency in the observer frame corresponding to 
electron LF 7 ^; the electron distribution function starts to fall off 
faster than 7 e~^ for 7 ' > 7 ^. The specific flux (flux per unit fre¬ 
quency) at Up due to PAS electrons is (e.g. Rybicki and Lightman, 
1979) 


fcs{Up) 


q^Bl,rNe{>i)l + Z 

meC^ ‘indj^ 


1 -h a 

meO’I'^R 4:ndj^ 


1 



7 ' 


p < 1 

(41) 

p > 1 


where Ne{> ip) is the total number of electrons inside PASs with 
LF Xi ip (which is ti times the integral of the source function 
given in equations|^& 37 1 , ti is the average time electrons spend 
in acceleration regions, and is the luminosity distance of the 
source at redshift 2 . 

The synchrotron spectrum due to electrons outside PASs is 
either /„ oc u~^^^ or depending on whether p 

is larger or smaller than 1 , and the ordering of u and synchrotron 
characteristic frequencies. 

The synchrotron flux at Up due to electrons outside PASs can 
be estimated using the distribution function calculated in the previ¬ 
ous subsection (eqs. |35[[^|38[|39y For the case we are considering 
where the guide field strength is of order B'q, the magnetic fields 
outside and inside PAS are of similar strength, and in that case the 
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fdziyp) ~ fcsiyp) 


t coot ifp) 

t'cs 


fcsfp) 


t'dzfo 

tcs^p 


(42) 


where fcs{vp) is synchrotron flux at Vp due to electrons inside 
PASs (see eq. [4T] (, and 


tcool (Tp) 


flTrrTleC 

(^TB'j-i'p 


{2.Zs)L^jfp-fTlRl^ 


(43) 


is synchrotron cooling time for an electron of LF 7 ^ outside PASs. 
Therefore, the ratio of synchrotron flux at Vp due to electrons inside 
and outside PASs is 


_ fcsjl'p) 

“ Uz{vp) 


^cool 


iJ'p)' 


(44) 


At a frequency v between Vmin. and Vp (assuming that idmin > r'c), 
the ratio of the flux from the two regions i^ 


fcsf) 

fdz{y) 


TZp 




l/3<p< 1 

p > 1 


(45) 


And the flux ratio at a frequency such that Vo <v < Vmin is 


fcsf) 

fdz{v) 


TZp 



(2-p)/2 


1/2 


1/3 <p < 1 

(46) 

p > 1 


The observed spectrum is a superposition of synchrotron ra¬ 
diation from electrons in PASs and electrons outside acceleration 
regions. We have provided all the relevant equations to determine 
specific flux at an arbitrary frequency from these two sources. We 
note that if the observed flux between Vmin and Vp is dominated 
by synchrotron radiation within PASs then the spectrum would be 
fi, oc which is harder by than the case where the 

flux comes mostly from electrons in the medium between PASs. 
Figure]^ shows the relative contributions of the two sources (upper 
panel), and the spectral index of the observed flux (lower panel), 
as a function of frequency. The figure clearly shows that the flux at 
Vp is dominated by electrons in PASs and thus the spectral index is 
harder. Electrons outside PASs make the dominant contribution to 
the observed flux at sufficiently low frequencies (below VpjXtP for 
the parameters considered in Fig. |^, and therefore the spectrum 
is softer. This behavior — softening of spectrum with decreasing 
frequency — is opposite to what we observe in the synchrotron ra¬ 
diation from particles that are accelerated in shocks. This spectral 
feature may be a way to determine if magnetic or shock dissipation 
of jet energy is responsible for the observed radiation. 

In practice, the spectral softening before the peak could be 
produced in a baryonic jet if there are multiple emission processes 
at work, for instance synchrotron radiation with a thermal compo¬ 
nent or synchrotron self-Compton (SSC). For a synchrotron - 1 - ther¬ 
mal component, the expected softening at low energies is very large 
which cannot be confused with the softening we expect for recon¬ 
nection described above. It is more difficult to tell apart a kinetic 
jet with SSC radiation mechanism and a magnetic jet where syn¬ 
chrotron dominates. The SSC should produce two distinct peaks 


° and Vo are synchrotron frequencies in the observer frame for elec¬ 

trons of Lorentz factors 7 ' . and 7 ' respectively in a magnetic field of 
strength Bq ; 7 / is given by eq. |30| and 7 /^^^ is the electron LF below 
which the average distribution function inside PASs either drops off or rises 
less rapidly than 7 /”^. 




Figure 4. The upper panel shows the ratio of synchrotron flux from elec¬ 
trons in acceleration regions (PASs) and electrons in regions outside; three 
different lines correspond to three different values of p: 0.5 (solid line), 
1.0 (dotted line) and 1.5 (dashed line). The parameters for these calcula¬ 
tions are: L = 10^®erg/s, F = 10 , cr = 10^, R = lO^^cm, eq = 0.1, 

= 0.1,7^ = 'y'rr.axn'min = Tmaa: /10^. and = 10-®XB/(cr). 

The lower panel shows the spectral index for the observed flux, i.e. a = 
dln{fcB + fdz)l<ttav for the same three values of p. Note that for the 
parameters chosen for these calculations the obseived flux is dominated by 
electrons radiating in PASs for v -J 10~^Vp and therefore the spectral index 
in this frequency range is a ~ — (p — l)/ 2 , i.e. the spectrum is hard. 


in the spectrum, one for synchrotron and one for inverse Comp¬ 
ton, and in this case the radiation at lower energies dominated by 
the synchrotron process can be softer. The spectrum for the mag¬ 
netic model on the other hand has just one peak, and therefore in 
principle it can be distinguished from the SSC model. However, in 
practice, the first peak for the SSC model may be at a frequency 
that is below the observing band, and that would make the task 
of identifying a Poynting jet more difficult. There are examples of 
low energy spectral softening in astrophysical objects, e.g., GRB 
090926A (Ackermann et al. 2011), but it is difficult to say whether 
or not it is due to multiple components to the spectrum or from the 
spectral feature described in this paper. We looked for the spectral 
feature we predict for a magnetic reconnection model in solar flares 
where magnetic dissipation is widely believed be at work (e.g. Lin 
et. al, 2003). Unfortunately, solar flares have a large thermal com¬ 
ponent, so the presence of a soft spectral feature at low frequencies 
is difficult to discern in these transients. 

The ratio of the luminosity in synchrotron and IC radiations 
is equal to the ratio of energy densities in magnetic field and pho¬ 
tons. Since only a small fraction of energy in magnetic fields in 
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dissipated in current sheetQ it is expected that the IC luminosity 
of a Poynting jet would be smaller than the synchrotron luminosity. 
The IC spectrum, or to be precise synchrotron-self-Compton spec¬ 
trum, is straightforward to calculate using the results for particle 
distribution and synchrotron spectrum described above. 


2.4 Constraint on the number of current sheets in the jet 


The radius interval where conversion of magnetic energy to ther¬ 
mal energy for a Poynting jet takes place depends on the magnetic 
field configuration and instabilities that develop in the jet. These 
are very difficult to calculate with any confidence. However, some 
general considerations described in this subsection provide broad 
guidance which can be used to constrain the dissipation radius and 
the number of current sheets in the jet. 

Consider a short segment of the jet that was launched at ra¬ 
dius f?o- The dissipation of magnetic energy in this segment could 
take place anywhere between Rq and ~ Rd (the deceleration radius 
of the jet), either gradually over this long distance interval or sud¬ 
denly within a short distance. Once the reconnection gets started 
at one location in the jet, it could trigger magnetic field dissipa¬ 
tion at other sites — possibly as a result of plasma outflow from 
this region or magnetic field reconfiguration propagating at Alfven 
speed and triggering reconnection at other sites — and these sec¬ 
ondary reconnection sites lie in a region of the jet that is in causal 
contact with the current sheet triggering these events. It is unlikely 
that most of the energy of this segment of the jet under consider¬ 
ation will be dissipated at a radius smaller than {St)cT^ because 
of causality considerations (provided of course that different parts 
of this segment of the jet don’t independently develop instability 
and/or reconnection centers). We can describe the dissipation with 
radius as a monotonically increasing function of radius, (^b{R)',Cb 
is the fraction of the magnetic energy in the segment that is dissi¬ 
pated or converted to bulk kinetic energy of the jet. 

For reconnections in a jet consisting of stripped magnetic wind 
geometry (magnetic fields reversing direction over distance of ro 
in the lab frame), (^b oc R^^^ (e.g. Drenkhahn and Spruit, 2002; 
Kumar & Zhang, 2014), and the proces s is completed at a radius 
Rc ~ roT^/eo; where eo as defined in i 2.1 
and magnetic fields and Vp ^ 


is the ratio of electric 


eoc is the speed for plasma flow into 
current sheets. So, although, the magnetic field dissipation process 
in this case is slow and extends over a large distance interval of 
Ro — roF^/eo 3> Ro, roughly half of the magnetic energy is in 
fact dissipated within a factor of a few of roP^/eo. This follows 
from causality, i.e. the size of the region where magnetic field is 
dissipated cannot increase at a speed faster than light, and hence 
the radial width of region where field has been dissipated grows 
proportional to R/F (in jet comoving frame) and that is the reason 
that a good fraction of magnetic energy dissipation occurs within a 
factor a few of the terminal radius where the dissipation process is 
completed. This property is likely to be generic, and independent 
of magnetic field geometry and reconnection model. 

Current sheets are likely to form and disappear on a time 
short compared with the dynamical time (R/ (cT) in jet comov¬ 
ing frame). We envision that there are an average Hs PASs present 
in the region at any given time, for a time duration ~ R/ (cT), and 


^ In general, it is highly unlikely for magnetic fields on the opposite sides 
of currents sheets to be exactly anti-parallel, and that limits the efficiency 
for converting magnetic energy to particle energy and radiation. 


the average length of these current sheets in jet comoving frame is 

a 

The total volume of plasma in the jet in the causally connected 
region at R is 

V' ~ (47) 

provided that the jet opening angle is > Qj. Since plasma flows 
into current sheets at speed v'p ~ cqC, the total volume of plasma 
passing through current sheets in a dynamical time is 

Vplasma,cs ~ CO (R/F) . (48) 

If the fraction of the magnetic energy in the jet dissipated in this 
region is ^b{R), then that means that the total volume of plasma 
passing through current sheets in volume 12' should be ()bV/. Thus, 
we obtain the number of PASs in the region to be 

~ [CB(/?)/eo]{i?/(r£;)}". (49) 

A lower limit for Ks can be obtained by substituting ~ R/F 
in equation l |49| l, which gives 

H.>Cs(i?)/£o. (50) 


And a generous upper limit on the number of current sheets can 
be obtained by taking £( ~ £a (the distance an electron trav¬ 
els in order to get accelerated to LF Xax)- Using i27 i we find 
5 X This upper limit is much too 

large to be of practical use. The length of an acceleration region 
can be much larger than £a when electron acceleration is balanced 
by radiative losses, and in that case far fewer number of PASs 
are needed to process magnetic energy to radiation. Let us take 
£'s = rjl'a, with rj ~ 10 ® that is needed to ensure that the observed 
radiation is dominated by electrons in PASs (as opposed to elec¬ 
trons in inter-PAS regions) and therefore the emergent spectrum is 
hard (see 92.3|. This results in 




-5x 


(51) 


3 DISCUSSION 

This work was motivated in part by a puzzle regarding gamma-ray 
bursts. A broad class of models for 7 -ray emission from GRBs is 
based on the jet being baryonic, which moves with a Lorentz fac¬ 
tor > 10^. The kinetic energy of baryons in the jet is converted to 
particle thermal energy via a series of shocks, and radiated away 
via the synchrotron process. According to this model, the spectrum 
below the peak should be oc or softer whereas the ob¬ 

served spectra for most bursts are close to i.e. much harder than 
the baryonic jet model predicts (e.g. Ghisellini et ah, 2000; Kumar 
& McMahon, 2008). The origin of this problem can be traced to 
the fact that particles are accelerated while crossing the shock front 
but otherwise they cool rapidly as they travel down-stream. There¬ 
fore, the number of electrons increases rapidly with decreasing LF 
(dn'f,/dj'f, oc 7 e~^ or faster) and that is the reason for the soft spec¬ 
trum for a generic model that is based on dissipation of baryonic 
jet energy in shocks. 

What we find is that if the GRB jet were not baryonic but 
Poynting, then the dissipation of magnetic fields and particle accel¬ 
eration provides a way out this problem. This is because particles 
can be kept in acceleration regions for a time much longer than the 
their radiative cooling time, thereby preventing the development of 
a large population of lower energy electrons that give rise to a soft 
spectrum. The spectrum of radiation from electrons in the region 
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in between PASs (where electrons are undergoing cooling without 
acceleration) is soft like that it is for the shock model, but the spec¬ 
trum emanating from PASs is hard because the powerlaw index for 
the particle distribution function in current sheets has p « 1. We 
have shown in §2.3| that the observed spectrum, which is a superpo¬ 
sition of contributions from the two regions (PASs and inter-PASs), 
is hard when the average time electrons spend in acceleration re¬ 
gion is much larger than their synchrotron cooling time. 

One of the general results reported here may be able to deter¬ 
mine whether a jet is baryonic or Poynting — for a Poynting jet, the 
spectrum below the peak softens with decreasing frequency, which 
is opposite to the case of a baryonic jet where shocks convert jet 
energy to radiation via the synchrotron process. 
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